Spatial normalization for voxel-based lesion symptom mapping: impact of registration approaches

BackgroundVoxel-based lesion symptom mapping (VLSM) assesses the relation of lesion location at a voxel level with a specific clinical or functional outcome measure at a population level. Spatial normalization, that is, mapping the patient images into an atlas coordinate system, is an essential pre-processing step of VLSM. However, no consensus exists on the optimal registration approach to compute the transformation nor are downstream effects on VLSM statistics explored. In this work, we evaluate four registration approaches commonly used in VLSM pipelines: affine (AR), nonlinear (NLR), nonlinear with cost function masking (CFM), and enantiomorphic registration (ENR). The evaluation is based on a standard VLSM scenario: the analysis of statistical relations of brain voxels and regions in imaging data acquired early after stroke onset with follow-up modified Rankin Scale (mRS) values.Materials and methodsFluid-attenuated inversion recovery (FLAIR) MRI data from 122 acute ischemic stroke patients acquired between 2 and 3 days after stroke onset and corresponding lesion segmentations, and 30 days mRS values from a European multicenter stroke imaging study (I-KNOW) were available and used in this study. The relation of the voxel location with follow-up mRS was assessed by uni- as well as multi-variate statistical testing based on the lesion segmentations registered using the four different methods (AR, NLR, CFM, ENR; implementation based on the ANTs toolkit).ResultsThe brain areas evaluated as important for follow-up mRS were largely consistent across the registration approaches. However, NLR, CFM, and ENR led to distortions in the patient images after the corresponding nonlinear transformations were applied. In addition, local structures (for instance the lateral ventricles) and adjacent brain areas remained insufficiently aligned with corresponding atlas structures even after nonlinear registration.ConclusionsFor VLSM study designs and imaging data similar to the present work, an additional benefit of nonlinear registration variants for spatial normalization seems questionable. Related distortions in the normalized images lead to uncertainties in the VLSM analyses and may offset the theoretical benefits of nonlinear registration

Bringing Anatomical Information into Neuronal Network Models

For constructing neuronal network models computational neuroscientists have access to wide-ranging anatomical data that nevertheless tend to cover only a fraction of the parameters to be determined. Finding and interpreting the most relevant data, estimating missing values, and combining the data and estimates from various sources into a coherent whole is a daunting task. With this chapter we aim to provide guidance to modelers by describing the main types of anatomical data that may be useful for informing neuronal network models. We further discuss aspects of the underlying experimental techniques relevant to the interpretation of the data, list particularly comprehensive data sets, and describe methods for filling in the gaps in the experimental data. Such methods of `predictive connectomics' estimate connectivity where the data are lacking based on statistical relationships with known quantities. It is instructive, and in certain cases necessary, to use organizational principles that link the plethora of data within a unifying framework where regularities of brain structure can be exploited to inform computational models. In addition, we touch upon the most prominent features of brain organization that are likely to influence predicted neuronal network dynamics, with a focus on the mammalian cerebral cortex. Given the still existing need for modelers to navigate a complex data landscape full of holes and stumbling blocks, it is vital that the field of neuroanatomy is moving toward increasingly systematic data collection, representation, and publication

26th Annual Computational Neuroscience Meeting (CNS*2017): Part 3 - Meeting Abstracts - Antwerp, Belgium. 15–20 July 2017

This work was produced as part of the activities of FAPESP Research,\ud Disseminations and Innovation Center for Neuromathematics (grant\ud 2013/07699-0, S. Paulo Research Foundation). NLK is supported by a\ud FAPESP postdoctoral fellowship (grant 2016/03855-5). ACR is partially\ud supported by a CNPq fellowship (grant 306251/2014-0)

25th annual computational neuroscience meeting: CNS-2016

The same neuron may play different functional roles in the neural circuits to which it belongs. For example, neurons in the Tritonia pedal ganglia may participate in variable phases of the swim motor rhythms [1]. While such neuronal functional variability is likely to play a major role the delivery of the functionality of neural systems, it is difficult to study it in most nervous systems. We work on the pyloric rhythm network of the crustacean stomatogastric ganglion (STG) [2]. Typically network models of the STG treat neurons of the same functional type as a single model neuron (e.g. PD neurons), assuming the same conductance parameters for these neurons and implying their synchronous firing [3, 4]. However, simultaneous recording of PD neurons shows differences between the timings of spikes of these neurons. This may indicate functional variability of these neurons. Here we modelled separately the two PD neurons of the STG in a multi-neuron model of the pyloric network. Our neuron models comply with known correlations between conductance parameters of ionic currents. Our results reproduce the experimental finding of increasing spike time distance between spikes originating from the two model PD neurons during their synchronised burst phase. The PD neuron with the larger calcium conductance generates its spikes before the other PD neuron. Larger potassium conductance values in the follower neuron imply longer delays between spikes, see Fig. 17.Neuromodulators change the conductance parameters of neurons and maintain the ratios of these parameters [5]. Our results show that such changes may shift the individual contribution of two PD neurons to the PD-phase of the pyloric rhythm altering their functionality within this rhythm. Our work paves the way towards an accessible experimental and computational framework for the analysis of the mechanisms and impact of functional variability of neurons within the neural circuits to which they belong

Mapping causal functional contributions derived from the clinical assessment of brain damage after stroke

Lesion analysis reveals causal contributions of brain regions to mental functions, aiding the understanding of normal brain function as well as rehabilitation of brain-damaged patients. We applied a novel lesion inference technique based on game theory, Multi-perturbation Shapley value Analysis (MSA), to a large clinical lesion dataset. We used MSA to analyze the lesion patterns of 148 acute stroke patients together with their neurological deficits, as assessed by the National Institutes of Health Stroke Scale (NIHSS). The results revealed regional functional contributions to essential behavioral and cognitive functions as reflected in the NIHSS, particularly by subcortical structures. There were also side specific differences of functional contributions between the right and left hemispheric brain regions which may reflect the dominance of the left hemispheric syndrome aphasia in the NIHSS. Comparison of MSA to established lesion inference methods demonstrated the feasibility of the approach for analyzing clinical data and indicated its capability for objectively inferring functional contributions from multiple injured, potentially interacting sites, at the cost of having to predict the outcome of unknown lesion configurations. The analysis of regional functional contributions to neurological symptoms measured by the NIHSS contributes to the interpretation of this widely used standardized stroke scale in clinical practice as well as clinical trials and provides a first approximation of a ‘map of stroke’

Reaction-diffusion models in weighted and directed connectomes.

Connectomes represent comprehensive descriptions of neural connections in a nervous system to better understand and model central brain function and peripheral processing of afferent and efferent neural signals. Connectomes can be considered as a distinctive and necessary structural component alongside glial, vascular, neurochemical, and metabolic networks of the nervous systems of higher organisms that are required for the control of body functions and interaction with the environment. They are carriers of functional phenomena such as planning behavior and cognition, which are based on the processing of highly dynamic neural signaling patterns. In this study, we examine more detailed connectomes with edge weighting and orientation properties, in which reciprocal neuronal connections are also considered. Diffusion processes are a further necessary condition for generating dynamic bioelectric patterns in connectomes. Based on our precise connectome data, we investigate different diffusion-reaction models to study the propagation of dynamic concentration patterns in control and lesioned connectomes. Therefore, differential equations for modeling diffusion were combined with well-known reaction terms to allow the use of connection weights, connectivity orientation and spatial distances. Three reaction-diffusion systems Gray-Scott, Gierer-Meinhardt and Mimura-Murray were investigated. For this purpose, implicit solvers were implemented in a numerically stable reaction-diffusion system within the framework of neuroVIISAS. The implemented reaction-diffusion systems were applied to a subconnectome which shapes the mechanosensitive pathway that is strongly affected in the multiple sclerosis demyelination disease. It was found that demyelination modeling by connectivity weight modulation changes the oscillations of the target region, i.e. the primary somatosensory cortex, of the mechanosensitive pathway. In conclusion, a new application of reaction-diffusion systems to weighted and directed connectomes has been realized. Because the implementation was realized in the neuroVIISAS framework many possibilities for the study of dynamic reaction-diffusion processes in empirical connectomes as well as specific randomized network models are available now

A spiking network model explains multi-scale properties of cortical dynamics

Neural networks in visual cortex are structured into areas, layers, and neuronal populations with specific connectivity at each level. Cortical dynamics can similarly be characterized on different scales, from single-cell spiking statistics to the structured patterns of interactions between areas. A challenge of computational neuroscience is to investigate the relation of the structure of cortex to its dynamics. Network models are promising tools, but for technical and methodological reasons, they have been restricted to detailed models of one or two areas or large-scale models that reduce the internal structure of areas to a small number of differential equations.We here present a multi-scale spiking network model of all vision-related areas of macaque cortex that represents each area by a full-scale microcircuit with area-specific architecture based on a model of early sensory cortex [1]. The layer- and population-resolved network connectivity integrates axonal tracing data from the CoCoMac database with recent quantitative tracing data, and is systematically refined using dynamical constraints [2]. Gaps in the data are bridged by exploiting regularities of cortical structure such as the exponential decay of connection densities with inter-areal distance and a fit of laminar patterns versus logarithmized ratios of neuron densities. Simulations reveal a stable asynchronous irregular ground state with heterogeneous activity across areas, layers, and populations. In the presence of large-scale interactions, the model reproduces longer intrinsic time scales in higher compared to early visual areas, similar to experimental findings [3]. Activity propagates preferentially in the feedback direction, mimicking experimental results associated with visual imagery [4]. Cortico-cortical interaction patterns agree well with fMRI resting-state functional connectivity [5]. The model bridges the gap between local and large-scale accounts of cortex, and clarifies how the detailed connectivity of cortex shapes its dynamics on multiple scales